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' Recent developments in the study of shape-invariant Hamiltonians are briefly summa- 

' rized. Relations between certain exactly solvable problems in many-body physics and 

^ shape-invariance are explored. Connection between Gaudin algebras and supersymmet- 

, ric quantum mechanics is pointed out. 

o 

I 1 Introduction 

Supersymmetric Quantum Mechanics (SSM) is the name given to the study of 
particular pairs of Hamiltonians SSM can be motivated by considering the 
ground ground state wavefunction, tpoix), for a one-dimensional bound system. 
. Since tpoix) has no nodes it can be written as 



O : Mx) = exp ( / Wix)dx ] , (1) 

ff^ i where the function W{x) is related to the potential energy of the system. Introduc- 

' ing the operators 

A = W{x) + —^p, 
V2m 

^ : = W(x) - -l=p, (2) 

pH ^ V2m 

one can write the Hamiltonian of the system as 

X ■ H-Eo= iU, (3) 

, where Eq is the ground state energy. The ground state wavefunction satisfies the 

condition 

i|Vo) = 0. (4) 
It is straightforward to show that the supersymmetric partner potentials 

Hi = 

H2 = AA^ (5) 

have the same energy spectra except the ground state of Hi , the energy of which is 
zero. Potentials corresponding to these Hamiltonians are 

viix) = mxr ' 



/2m dx 

V2{x) = [Wix)r + ^^. (6) 
V2m dx 
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The partner potentials in Eq. (6) are called shape-invariant^ if they can be obtained 
from one another by changing their parameters: 

V2{x;ai) = Vi{x;a2)+R{ai), (7) 

where 02 is a function of ai, and the remainder R{ai) is independent of x. Eq. (7) 
is equivalent to the operator relation 

i(ai)it(ai) = it(a2)i(a2) + R{ai). (8) 

1.1 Algebraic Approach 

Shape- invariance problem was formulated in algebraic terms in Ref. [4]. In this 
formulation one introduces an operator which transforms the parameters of the 
potential: 

f{ai)0{ai)f-\ai)=0{a2). (9) 

Defining the operators 

B_=Bl = f'i{ai)A{ai) (10) 

one can show that the Hamiltonian can be written as 

H - Eo= A''A = B+B_. (11) 

Using the definitions given in Eq. (10), the shape-invariance condition of Eq. (8) 
takes the form 

[B_,B+]=R{ao), (12) 

where R{ao) is defined via 

R{an) = f (ai)i?(a„_i)f t(ai). (13) 
In terms of these new operators Eq. (4) takes the form 

B_|Vo) = 0, (14) 

i.e. the ground state is annihilated by the lowering operator 
One can easily establish the commutation relations^ 

[H, B^] = (i?(ai) + i?(a2) + • • +i?(a„))B!^ (15) 

[H, B!!] = -B!! (i?(ai) + i?(02) + • • +R{an)) . (16) 

i.e., B" IV-'o) is an cigcnstate of the Hamiltonian with the eigenvalue R{ai) + R{a2) + 
■ ■ +R{an)- The normalized eigenstate is 

IV'n) = ^ B+ ■ ■ ^ ^ =B+ J_^ E+|V>o). (17) 

^R{ai) + • • +i?(a„) + ^R{a^)+R{a2) WWi) 
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To identify the algebra we consider the commutation relations 

=i?(ao) (18) 

[B+,Riao)] = (i?(ai) - i?(ao))S+, (19) 

[B+, (i?(ai) - R{aa))B+] = {(i?(a2) - R{ai)) - (i?(ai) - R.{ao))}B+, (20) 



and so on. In general there are an infinite number of such commutation relations. If 
the quantities R{an) satisfy certain relations one of the commutators in this series 
may vanish. For such a situation the commutation relations obtained up to that 
point plus their complex conjugates form a Lie algebra with a finite number of 
elements. For example if the condition 

(i?(a2) - R(ai)) - {R{ai) - R{aa)) = (21) 

is satisfied then the algebra is** either SU{2) or SU {1,1). Most of the exactly 
solvable one-dimensional problems in quantum mechanics can be described by this 
algebra ^ . It can be shown that this algebra also describes for example both the 
bound and scattering states of the Poschl- Teller potential ^ as well as associated 
transfer matrices. 

1.2 Outlook on future applications 

Almost all exactly solvable one-dimensional potential problems encountered in quan- 
tum mechanics textbooks are shape invariant where the parameters are related by 
a translation ^ 

a2=ai+r]. (22) 

It should be emphasized that shape-invariance is not the most general integrability 
condition one can write for such potentials as there are exactly solvable problems 
which are not shape invariant^. There is a second class of shape invariant potentials 
where the parameters of the partner potentials are related by a scaling 

a2 = qai. (23) 

In this latter class, corresponding one-dimensional potentials are defined implicitly, 

but explicit expressions are not given. 

In searching for integrable models in two-dimensional statistical mechanics a 
relationship was uncovered between those models, three-dimensional Chern-Simons 
gauge theory and quantum groups These models, being completely integrable, 
can be written in a shape- invariant way corresponding to a shift in the parameters 

02 = qai + r). (24) 

The associated algebras arc called up-down algebras These developments suggest 
that there may be shape-invariant potentials where the parameters are related by 
linear- fractional transformations : 

a2 = {qai + r?)/(ai + rj') (25) 
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This is a completely unexplored direction of research as nothing is known about 

such intcgrablc systems. Recall that the notation ai, a2, etc. may represent not 
only single parameters, but also a set of them. In general one may suggest to simply 
relate these parameters by the transformation 

f{ai)0{ai)f-\ai)=0{a2). (26) 

where T is an element of any group, not just of SL(2,R) as suggested by the linear- 
fractional transformation and its limits that were so far employed. What kind of 
exactly solvable problems do wc obtain? At the moment this is an open question. 

The basic philosophy of this approach is to consider the parameters of the 
Hamiltonians as auxiliary dynamical variables. This is reminiscent of the path 
leading to the Interacting Boson Model To describe the quadrupole collectivity 
in nuclei one needs to consider a five-dimensional space. It is possible to formulate 
this problem in terms of boson variables however the problem is nonlinear written 
in terms of quadrupole bosons. By considering a parameter of the problem (boson 
number) as an additional degree of freedom. Interacting Boson Model introduced a 
scalar boson as a dynamical variable. This has led to the subsequent realization^^ 
of s and d bosons as pairs of nucleons coupled to the angular momentum L = and 
L = 2 

So far we talked about considering parameters of the shape-invariant problem 
as auxiliary dynamical variables. One can imagine an alternative approach of clas- 
sifying some of the dynamical variables as "parameters" . An example of this is 
provided by the supersymmetric approach to the spherical Nilsson model of single 
particle states The Nilsson Hamiltonian is given by 

iJ = ^ a|ai - 2fcL.S + fci/L^ (27) 

i 

Introducing the variable 

Ft=^<7,aJ (28) 

i 

one can show that the "Hamiltonians" 

Hi = F^F = ajai - a.L (29) 

i 

and 

H2 = FF^ = J2 aia] + a.L (30) 

i 

can be considered as supersymmetric partners of each other The shape-invariance 
condition of Eq. (8) can be written as 

FF^f = F^'F + R, (31) 

where the remainder is 

i? = C7.L - 3/4, (32) 
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i.e. in this example the radial variables are considered as the main dynamical 
variables and the angular variables are considered as the "parameters" . 

A number of applications of shape-invariance are available in the literature. 
These include i) Quantum tunneling through supersymmetric shape-invariant po- 
tentials^^; ii) Study of neutrino propagation through shape- invariant electron den- 
sities^^; iii) Investigation of coherent states for shape-invariant potentials ^^'^''; and 
iv) As attempts to devise exactly solvable coupled-channel problems, generalization 
of Jaynes-Cummings type Hamiltonians to shape- invariant systems In this 

article we focus on the applications to many-body systems. 

2 Many-Body Hamiltonians 

One can ask if these methods can be used to search for exactly-solvable many-body 
systems. It has been shown that the concept of supersymmetric shape-invariance 
can be utilized to derive the energy spectrum of Calogero-Sutherland model Here 
we discuss an alternative approach and first write down multiple commutators for 
a shape-invariant Hamiltonian 

[H,B+]=R{ai)B+ (33) 

[[H, B+],B+] = (i?(ai) - R{a2))Bl (34) 
[[[H, B+],B+],B+] = (Riai) - 2R{a2) + R{a3))Bl (35) 
[[[[H, B+\,B+IB+IB+] = iR{ai) - 3i?(a2) + 3R{a3) - R{,ai))BX (36) 

and so on. We wish to address the possibility of defining an exactly solvable problem 
through these commutation relations. We will consider B+ as a raising operator. 
We assume that the Hamiltonian H may or may not be in the form given by Eq. 
(11). We consider a generalized pairing problem with 

B+ = ^c,-5;. (37) 

3 

In Eq. (37) the pair creation operator in a single-j shell is defined as 

m 

where aj^ is the particle creation operator. If we assume that the shape-invariant 
Hamiltonian has only one- and two-body terms the commutator [[/f, -B+], B+\ will 
only involve products of four creation operators. Consequently the next nested 
commutator will vanish: 

[[[H,B+],B+],B+]=Q (39) 

Higher nested commutators will also vanish. This will place strong constraints on 
i?(a„), i.e. 

i?(a3) = -R{ai) + 2R{a2), (40) 
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R{a4) = R{ai) - 3i?(o2) + 3R{a3) (41) 

and so on. Consequently we can immediately write the energy eigenvalues and 
eigenstates of the Hamiltonian 

HBlliJo) = (nRiai) + ^Wn{n - 1)) |Vo), (42) 

where 

W = Ria2)-R{ai). (43) 
A similar approach was first given by Talmi . 

3 Connection to Gaudin Algebras 

The pairing model with a constant two-body interaction was solved exactly by 

Richardson . In a parallel development Gaudin developed an algebraic approach 
to solve many-body spin Hamiltonians Here we will explore the relationship 

between Gaudin's methods, algebraic methods developed to search for quasi-exactly 
solvable models and supersymmetric quantum mechanics. 

Following the notation of Ref. [29] we consider the function defined as 

N 

*(A)=n(A-Ci)e--/''^'\ (44) 

i 

where W{\) is an arbitrary function of A and are numbers to be determined. 
Introducing the operators 

A = W + ip, = W-ip, (45) 
it can be shown that the function defined in Eq. (44) satisfies the equation 



oV 1 gV W{\) 

^(A-Ci)(A-^,) 4"(A-^,] 



^. (46) 



Requiring the residue at to vanish yields the Bethe-ansatz conditions: 

Wm = Y.-^. (47) 

Inserting Eq. (47) into Eq. (46) we obtain 



(48) 



Provided that their superpotentials satisfy the condition given in Eq. (47), factor- 
ized supersymmetric Hamiltonians satisfy Eq. (48). Note that the right side of Eq. 
(48) in general depends on A, hence we cannot interpret the term that multiplies 
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the function ^ as an energy eigenvalue. However, for a number of limited cases 

(certain functions W(A) such as those that correspond to a harmonic oscillator) 
this A dependence drops out and one can recover the standard expressions for the 
energy eigenvalues^'^. 

The three generators of Gaudin's algebra (Jo (A), J±{X)) can be defined through 
the commutation relations 

[MX),JM = -^^^^^^^, (49) 

[J_(A),J4m)] = -2^^^4^^, (50) 

and 

[J±,o(A),J±,o(/i)] =0, (51) 

where A is, in general, a continuous parameter. Gaudin studied the eigenstates of 
the "Hamiltonian" ^'^ 

H{X) = Jo(A)Jo(A) - ij_(A)J+(A) - ij+(A)J_(A) (52) 

If a state | 0) which is annihilated by all J_ (A) can be identified 

J_(A)|0)=0, (53) 
then W{X) is introduced as the eigenvalue of Jo(A) on that state: 

Jo (A) I 0) = W^(A) I 0). (54) 
One can then find the eigenvalues and eigenstates of the "Hamiltonian" of Eq. (52): 

iJ(A) I <P) = E{X) I (55) 

where the eigenstates are 

I $) = J+(^;v) J+(Cjv-i) • • • J+(a) I 0), (56) 
and the eigenvalues are 

E{X) = W^X) + ly'(A) + 2^ (^ ^^^] ~ J^^'^ ^ . (57) 
In deriving the above equations the conditions 

W^te) = -E^^' ^'J = l'---'^ (58) 

were assumed to be fulfilled. 

The strategy of using Richardson-Gaudin methods to deal with many-body 
problems were employed by a number of authors Clearly there is a map- 
ping between the solutions of the Gaudin problem (Eq. (55)) and those of the 



factorized supersymmetric Hamiltonians (Eq. (48)). One may ask if this correspon- 
dence can be exploited to study pairing and related problems. 

The pair creation operator^® in a single-j shell is defined in Eq. (38), 



(59) 



its Hermitian conjugate, and the number operator span an SU (2) algebra (the so- 
called quasi-spin algebra). One can obtain a Gaudin algebra from the quasi-spin 
algebra by defining 

^+W = Ea^' (60) 



(and similar formulas for the other elements). This realization of the Gaudin algebra 
can be very useful in many-body systems. As a simple example we consider a system 
with s and p bosons and define three operators that satisfy Gaudin's commutation 
relations 



X — as A — Q!p 
B_(A) = [B+(A)]t, 



and 



X- as 



hp + 3/2' 
X — a„ 



It is easy to show that as as ap the quantity i?+(A)i3_(A) reduces to 

1 



X — an 



(61) 
(62) 

(63) 
(64) 



where P4 is the 0(4) pairing operator. One can then study a Gaudin- type Hamil- 
tonian which generalizes this operator 



H{X) = Bo{X)Bo{X) - ^B_{X)B+{X)B+{X)B_{X). 



(65) 



Following steps above one can show that this Hamiltonian is associated with the 
one-dimensional potential 



V{x) = - 



2 \x — as 



X- a. 



(66) 



Similar ideas could conceivably be useful in dealing with other many-body systems. 
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